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Abstract. The effect of fluctuations on the conductivity of ion channels is 
investigated. It is shown that modulation of the potential barrier at the selectivity site 
due to electrostatic amplification of charge fluctuations at the channel mouth exerts 
fS |' a leading-order effect on the channel conductivity. A Brownian dynamical model of 

ion motion in a channel is derived that takes into account both fluctuations at the 
channel mouth and vibrational modes of the wall. The charge fluctuations are modeled 
as a short noise flipping the height of the potential barrier. The wall fluctuations 
are introduced as a slow vibrational mode of protein motion that modulates ion 
qq , conductance both stochastically and periodically. The model is used to estimate the 

' contribution of the electrostatic amplification of charge fluctuations to the conductivity 

. of ion channels. 
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1. Introduction 

Ion transport through the channels in cellular membranes underlies the electrical 
signal transduction and processing by living organisms. Accordingly ion channels, 
being natural nanotubes, control a vast range of biological functions in health and 
disease. The understanding of their structure-properties relationship is the subject 
of intensive, ever-growing, fundamental and applied research in biology, physics, and 
nanotechnology [U [2] . A central problem in studies of ion permeation through biological 
membrane channels is to understand how channels can be both highly selective between 
alike ions and yet still conduct millions of ions per second [3] . Indeed, selectivity between 
ions of the same charge implies that there exists a deep potential well for conducting 
type ions at the selectivity site of the channel. On the other hand such channels can 
pass up to 10 8 ions per second pi] corresponding to almost free diffusion. 

Significant progress has been made towards an understanding of this problem over 
the last few decades. In particular, the molecular structure of the KcsA potassium 
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channel [5] that discriminates between Na + and K + was determined by crystallographic 
analysis. Furthermore, by detecting the size of the structural fluctuations [6] and 
conformational changes [7], it has become possible to provide the experimental 
information needed for molecular modelling of the dynamical features of the observed 
selectivity and gating [HIE]. In particular, the minimum radius of the selectivity filter 
in KcsA is ~0.85A, which is to be compared with 1.33A for the ionic radius of K + , 
suggesting that flexibility of the filter is coupled to ionic translocation [10] . It has 
therefore become apparent that fluctuations in the channel walls plays a fundamental 
role in maintaining high conductivity in highly selective channels (see also Elber [TT]). 

Another important source of modulation of the electrostatic potential identified in 
earlier research [T2| [T3] relates to the interaction of the ion in the channel with charge 
fluctuations in the bath solutions. The effect of current fluctuations and noise on the 
channel entrance rates and on the channel conductivity was also considered in [HI fl5] . 
It becomes clear that fluctuations of the electrostatic potential within ion channels 
induced by various sources may provide a key to the solution of the central problems 
of permeation and selectivity. Models of such fluctuations have thus become one of the 
central topics of research on the permeability of ion channels. It is important to note 
that dynamical models of ion motion in the channel can also provide a link between 
studies of the permeability of open channels and channels gating. Notwithstanding 
recent advances, theoretical modeling of the dynamical features of ion channels is still 
in its infancy. In particular, little is yet known about the relative importance of the 
different dynamical mechanisms and sources of fluctuations in the ion channels. 

In our earlier work we have started to develop a dynamical model [HI [T7J [TU 
[T9l l20l l2T] of ionic conductivity through open channels. It takes into account the 
coupling of ion motion to vibrations of the wall [HI [T7J and to charge fluctuations 
at the channel mouth [19], [20]. Our goal is to derive a self-consistent model that allows 
for analytical estimation of the potential barrier at the selectivity site and for the effects 
of fluctuations on the conductivity of the channels. The starting point of our approach is 
a self-consistent quasi-analytical solution of the Poisson and Nernst-Planck equations in 
the channel, and in the bulk [21], allowing for accurate estimation of the current- voltage 
characteristics of ion channels [22] (see also [23J, [24] [25], [26] [27]). The electrostatic 
channel potentials resulting from these estimates can be further used to estimate relative 
contribution to the channel conductivity from the different sources of fluctuations. 

In this paper we introduce a model of ion permeation that takes into account 
dynamical effect of the charge fluctuations through the resultant shot noise, and we 
demonstrate that the latter has a leading-order effect on the transition probabilities. We 
show that the charge fluctuations at the channel mouth can be modeled as a flipping of 
the electrostatic potential at the selectivity site, which fluctuates between two maximum 
values at a rate corresponding to the random arrivals of ions at the channel mouth. The 
developed theoretical framework will allow us in the future to include into the model 
the modulation of the potential at the selectivity site due to hydration effects inside the 
channel. 
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A model of 3D Brownian dynamics simulation of ions in the bulk and inside the 
channel is described in Sec. |2j Using results of the 3D simulations in the bulk we present 
in Sec. |3]a reduced model of an ion moving in the channel and interacting with the wall 
vibrational modes and with charge fluctuations at the channel entrance. The model uses 
the channel potential derived from a self-consistent solution of the Poisson equation and 
the flipping rates of the potential barrier obtained from simulations of Brownian motion 
of ions in the bulk. In Sec. HI we estimate analytically the mean first passage time of 
the channel. These estimates are based on the assumption that barrier-crossing and 
barrier-fluctuations are correlated for a general form of the potential barrier [28]. The 
mean first passage time is calculated as a weighted sum of the escape time r_ over the 
low barrier AE , and the escape time r + over the high barrier AEi. Conclusions are 
drawn and future directions of research are outlined in Sec. EJ 

2. Brownian dynamics simulations in the bulk and inside the channel 

We consider the following Brownian dynamical model of the ion permeation (see sketch 
in the Fig. [1] (left)). The system is made of three compartments of equal size. The middle 
block constitutes the protein through which there is a cylindrical hole approximating 
the open channel. To model the ion's coupling to the vibrational modes of the channel, 
we introduce a moving segment of the protein wall that is elastically bound to the wall. 
The moving segment is charged and represents the selectivity site. The dynamics of the 
ions in the bath and channel, and of the moving segment, are modeled using Brownian 
Dynamics (BD) simulations, see eqs. (J2J)- (JHI) coupled to the Poisson equation (CQ). 

-V-(5(f)V0(r-))=p(f), (1) 

rriiXi + rriijiXi = F Cji + F sr>i + F Hji + y/2m i 'y i k B T^ i {t), (2) 

mjXj + rrij'jjXj = F C j + F sr j + F H>j + y / 2m/)p^£,- (t) , (3) 

mi + m 7 x = F C + F ch + F iWji sin(/3) + ^2m-fk B T^(t) , (4) 

M5R + MT5R + K5R = F iwM cos(/3) + y/2MTk B Tu(t) (5) 

Here r^ = Xi — Xj, = \fij\, mi, Xi and qi, mcii and \j2m i 'yikBT ^{t) are the mass, 
position, charge, friction coefficient and the stochastic force of the ith ion. The distance 
between ions i and j is r^. In this work, for simplicity, we restrict the analysis to two 
types of ions: the index i will correspond to Na + , while index j corresponds to CI -1 . The 
motion of the charged residual of mass M at the selectivity site is characterized by the 
displacement 5R in the vertical direction from the equilibrium position R and an elastic 
force oc K5R. Note that in general values of the effective mass and diffusion constant of 
the ion moving within the channel may deviate from the corresponding values in bulk 
due to the nontrivial effect of hydration in the channel. Coupling between the motion 
of the ion in the channel and the normal mode of the wall oscillations is accounted for 
by the term Fi W corresponding to the Coulomb interaction between ion and charge at 
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Figure 1. Left: Sketch of the model. The middle block represents the protein, 
through which there is a cylindrical hole approximating the open ion channel. The 
moving charged segment of the protein wall, assumed elastically bound to the wall, is 
shown by the filled squares. Negative ions are shown by the filled circles, and positive 
ions by the open circles. Right: Sketch of an ion moving along the channel axis. The 
conducting ion is shown by red circle. The difference between ion coordinate x and 
location of the selectivity site x ss is z. The vertical displacement of the selectivity site 
is SR. The angle between a line connecting the ion with the selectivity site and the 
vertical axis is (3. 



the selectivity site. Keeping only the axial component of the force for the ion motion in 
the channel, and the radial component of F iw for the oscillating wall, we have 

= Q<ij z F . = Qij 6R 

lw ' 1 47iee d 2 of 1 ' 47iee d 2 d 

Here d = yj ((R + 5R) 2 + z 2 ), z = x — x ss where x is the coordinate of the ion in 
the channel along the channel axis, and x ss is the location of the selectivity site. An 
additional coupling of the ion motion to fluctuations of the channel wall is effected via 
modulation of the channel potential by the moving wall. This is taken into account 
through the term F ch on the rhs of eq. (jlj). The ions in the bulk (see eqs. (J2J, 0) 
do not feel the channel potential. Instead, their motion is governed by the Coulomb 
interaction Fc, by the short-range interaction F sr , and by hydration forces Fg (see [3]) 

r C,i 2-uj=l Aweeorfj n s ' F * r t ^j=l ry ' 

Fh„ = Ef^^oexp (^) sin ( 2 vrta _ 



where A = yjl + (^-) 2 and a = arctan(^-). 

In addition the effect of the surroundings is taken into account by an average 
frictional force with frictional coefficient and a stochastic force y/ 2m i r ) i kBT£ ti (t) . 
The addition of the pairwise repulsive 1/r 9 soft-core interaction potential ensures that 
ions of opposite charge, attracted by the inter-ion Coulomb force, do not collide and 
neutralize each other. Uq and R c are respectively the overall strength of the potential 
and the contact distance between ion pairs. The oscillating part is added to the potential 
and takes explicitly into account the internuclear separation for the two solvents, where 
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a w is the oscillation length, a e is the exponential drop parameter, and Rij is the origin 
of the hydration force which is shifted from R c by +0.2 A for like ions and by -0.2 A 
otherwise [29]. F c h is the dielectric force in the channel, obtained by solving Poisson's 
equation numerically using finite volume methods (FVM) [30]. We use the Langevin 
equation to model the collective motion of the atoms forming the channel protein charged 
ring located at the selectivity filter. In this way, our analysis is based on the assumption 
that the movement of structural domains of the channel protein may be described as the 
motions of independent, elastically bound Brownian particles [31]. We have included 
the damping term MT5R and the corresponding random force y/2MTkBTu(t), whose 
amplitude is related to the damping constant via the fluctuation-dissipation theorem. 
The function u(t) is a Gaussian white noise. Q is the total fixed charge on the flexible 
ring interacting with an ion of charge qj on the channel axis z. R is the channel radius, 
5R is a small variation of the channel radius and K is the elastic constant of the channel 
protein. In the particular case of the Gramicidin A (GA) channel, the value of the 
elastic constant is estimated by calculating the root-mean-square deviation (RMSD) of 
the backbone forming its central part, together with the single ion potential of a Na + 
ion as it traverses the channel. Using data from [32], we obtain an elastic constant of 
~ 1.6577 N/m for a maximum RMSD of 1 A. The GA channel molecular weight M 
is about AkDa = 6.64 x 10 -24 kg [33]. The diffusion constant of the protein in the 
membrane is between 10 -14 and 10~ 16 m 2 /s [34] . 

The parameters of the ion-ion interaction are presented in Table. [H 

Table 1. Parameters used in the calculation of the short range ion- ion interaction 
with hydration 



Ions 


U [k B T] 


R C [A) 


R[A] 


Oio [A] 


Oe [A] 


Na-Na 


0.5 


3.50 


3.7 






Na-Cl 


8.5 


2.76 


2.53 


2.76 


1.4 


Cl-Cl 


1.4 


5.22 


5.42 







Other parameters used in simulations are: 

Dielectric constants: 62 = 80, S\ = 2; 

Masses (in kg): m Na = 3.8 x 10" 26 , m cl = 5.9 x 10" 26 ; 

Diffusion coefficients (in mV 1 ): D Na = 1.33 x 10~ 9 , D C i = 2.03 x 10~ 9 , 

(Note that D is related to the friction coefficient via D = ^ M ^)\ 

Ion radii (in A): r Na = 0.95, r C i = 1.81; 

Temperature: T = 298 K. 

3. Reduced model for ion channel conduction 

To derive the reduced model we notice that Eqs. (DO)-© correspond to a many-body 
problem with widely-varying timescales, ranging from ps (ion fluctuations) to fis (wall 
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vibrations). We further assume that the channel is occupied most of the time by only 
one ion, and that the transition rate of ions through the channel is mainly determined by 
escape over the potential barrier at the selectivity site. Then the effect of the many-body 
ion dynamics in the bulk on the ion motion in the channel is twofold: (i) a delivery of 
the ions to the channel mouth and (ii) modulation of the channel potential by the charge 
fluctuations at the channel mouth. Under these physiologically plausible assumptions 
one can separate the ion motion in the channel from the many-body ion dynamics in 
the bulk. The resulting equations of ion dynamics in the channel coupled to the wall 
fluctuations can be written as follows 

dV(x,t) 



rwyx 



dx 



+ ^2mjk B T^{t), (6) 



M5R + MT5R + K5R = F iW)M cos(/3) + ^2MTk B Tis(t) . (7) 

Note that the reduced motion of the conducting ion is overdamped while the wall 
fluctuations are damped. Damped vibrational mode models the relatively slow (on 
a time scale of ns) motion of the protein of the channel walls that was suggested 
to be essential for the ion transport process. 




Figure 2. Simulations of Brownian dynamics in the bulk liquid. Left: The number 
of ions as a function of the distance between ions of different types: C=400mMol; 
Box size: 40 A. Forces included in the simulations were: (i) Coulomb interaction; 
(ii) short range repulsion; (iii) hydration. Right: The arrival time distribution for 
positives(solid line) and negative(dashed line) for cylindrical channel of radius r=6A 

The potential V(x,t) in our approximation has three main contributions: (i) the 
potential of Coulomb interaction with ions in the bulk solution Vc', (ii) the electrostatic 
potential induced by interaction with the channel protein V c h, (iii) the potential of 
Coulomb interaction with the wall charge at the selectivity site. By an averaging 
procedure, the effect of multi-ion motion in the bulk solutions is reduced to Coulomb 
interaction with ions arriving at the channel mouth. The later process can be viewed 
as a stochastic Poisson process or as shot noise that modulates the potential barrier for 
the conducting ion at the selectivity site. The goal of the present paper is to estimate 
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analytically the effect of this potential modulation on the channel conductivity as will 
be discussed in details in the next section. 

To quantify the effect of multi-ion motion in the bulk on the conducting ion at the 
selectivity site, we have simulated ion's Brownian dynamics in the bulk. The resulting 
ion-ion distributions in the bulk are shown in the Fig. [2] (left). We emphasize that these 
distributions are very close to those obtained earlier in both BD simulations [3] and 
experiments [36]. The arrival time distributions for Na + and Cl~ ions at the channel 
mouth (defined as a cylindrical section of radius R and length R) obtained in our 
simulations is shown in Fig. [2] (right). Note that these distributions are exponential for 
both Cl~ and Na + with mean arrival times r~ = 471 ps and r + = 365 ps respectively. 
These estimates are in agreement with the theoretical estimates in [37] 

1 

Tarrival ~ (8) 

where c is the ion concentration and D = is the diffusion coefficient. 

my 




Figure 3. (left) Charge fluctuations at the channel mouth. The positive charge is 
shown by the solid line. The negative charge is shown by the dashed line, (right) The 
potential energy profiles as a function of the position of the ion when: the first ion is 
fixed at the channel mouth (dashed line) and the second is moving along the channel 
axis. The solid line corresponds to the potential energy on a single ion moving on 
the channel axis, and the potential energy of the passive channel (dotted line). The 
vertical dashed lines show the channel entrance and exit. The height of the potential 
energy barrier seen by a single ion at the selectivity site as it moves from left hand 
to right hand of the channel is denoted AEq. In the presence of a second ion at the 
channel's left mouth this barrier is reduced to AEi. 



The time evolution of the charge in the channel mouth is shown in Fig. [3] (left). It 
can be seen that the charge at the channel mouth is a Poisson process with the three main 
states +le, 0, and — le. As a first approximation it is convenient to divide the states 
of the channel potential affected by the charge fluctuations into two effective states: 
(i) state of high conductivity, corresponding to +le, and (ii) state of low conductivity, 
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corresponding to or — le charge at the channel mouth. In this approximation the effect 
of three states of the potential is taken into account by asymmetry of the transition rates 
between the two effective states. The corresponding transition rates can be estimated 
as a 1 * 1 = 1/ (T±), where (T±) are mean residence time of two effective states, giving 
(ct + ) _1 = 0.22 ns and = 0.38 ns. The occupation probabilities of these two states 

are 0.36 and 0.64 respectively. The effect of the wall oscillations on the channel potential 
was estimated earlier [20] and for simplicity the wall will be assumed rigid in the rest of 
the paper. 

To estimate the effect of charge fluctuations on the value of the channel potential 
we solve the Poisson equation for various positions of the conducting ion along the 
channel axis in two cases: (i) when there are no other ions at the channel entrances; (ii) 
when there is one positive ion at the left entrance to the channel. The results of these 
calculations are shown in Figure, fright). The low-conductivity effective state of the 
channel is shown by the blue solid line and corresponds to the potential barrier « l&^T 
at the selectivity site. The high conducting state is shown by the red dashed line and 
corresponds to a potential barrier height ~ 0.2/c#T, i.e. practically no barrier state. 
It can be seen that the charge fluctuations are enhanced in channels of low dielectric 
constant, resulting in strong modulation of the potential barrier at the selectivity site [20] 
(see also [38J). 




Figure 4. Approximation of the fluctuating potential. 



It is therefore possible to build a simple model capable of coupling the motion of 
ions in the channel to the bath solution. The channel potential becomes 

y _i_ v V+-V 
V(x,t) = ^- + ^ T - X (t), *(f) = ±l, (9) 

AE A E 
V + = -{x-x ) + E+, V^ = °-{x-x ) + E , 

where x(t) is a Poisson random force with two transition rates a± between the states 
+1 and —1 . The charge fluctuations at the channel mouth thus result in flipping of the 
potential. Here AE = E^ — Eq and AEi = E^ — Eq , with AE > AE\ are respective 
barriers of the potential in two states, and x m = X\ — x . 
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A direct analogy can be made between the model described by Eq. (JSJ) and the 
model described by Ziircher [28] whose barrier fluctuation is controlled by a dichotomic 
noise of zero mean and exponential correlation. An approximation of the fluctuating 
potential is sketched in Figure. HI The similarity of the two problems suggests that there 
is some possibility of semi-analytical estimations of the effect of the charge fluctuations. 



4. Estimation of the mean channel crossing time 

We are interested in the transition of a particle initially trapped at the channel selectivity 
filter. This corresponds to the motion of the ion across the interval [xo,^]. The 
approximate potential in this interval is given by ([9]). We are interested at the 
unidirectional current, so there is no backward flow of ions. We assume that on average, 
the channel is always occupied by a single Na + ion. This arises from the fact that, when 
an ion is sitting in the middle of the channel, it is almost impossible for a second ion to 
enter the channel due to the height of the barrier at the left entrance of the channel, as 
can be seem from Figure. HI 

As a first approximation, we assume that the mean first passage time (MFPT) is 
only determined by escape. Therefore the MFPT for the channel is expected to be a 
function of the two times r + and r_ corresponding to the escape times from the potential 
minimum in two effective states of the potential. Our derivation follows very closely the 
earlier discussion by Ziircher [28] with the difference that, here, we are interested in the 
asymmetric case with two transitions rates. 

Assuming no back flow, the quantities t±(x) are defined by (see [39J): 

-ldV + dr+ k B Td 2 r + + . . 

j 1— H YT ~ a T + + a r_ = -p+, (10) 

my dx dx my dx z 

-ldV-dr„ k B Td 2 T_ , 

-, -j— + — - - a r_ + a + r+ = (11) 

my dx dx my dx 

The potential jumps between positive and negative values, with respective rates a + and 
a~. p± are the occupation probabilities of these states. We choose a reflecting boundary 
condition (BC) at the bottom of the barrier x = Xq and absorbing BC at the top of the 
barrier x — x±, 

dr+(x = Xn) , , „ 

- = 0, t ± (x = xi =0, (12) 

dx 

With r + and r_ specified, the mean exit time of the Brownian particle that is trapped 
at the selectivity filter x = Xq is given by: 

r = r+(x ) + r_(x ), (13) 

Following Ziircher [22], the calculations of t± is straightforward. We introduce 

^ AE ' + AC °, A^ A£ '- Afl °, (14) 
2 2 

and 

a(x) = a + r + (x) + a~r_(x), S(x) = a + r + (x) — a~r_(x). (15) 
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Figure 5. Left: MFPT as a function of a + . Right: Current as a function of a + . 



We find the coupled differential equations 
-DE da d 2 a . , . , , Nr DA d5 

+ ^~r~ ; ) + \ a P+ + a P-) = \ a P+ ~~ a P-)° r~^~r-> I 1 ") 



x m kBT dx dx 2 ' XmksT dx' 

-DE d5 d 2 5 . , _ . . _ , DA da 



x m hj^T dx dx 2 x m kj^T dx 

where D = D Na and the boundary conditions, 

daix = Xo) . = Xn) 

<r(a; = xi) = 0, — —, = 0, 5{x = x 1 )=0, —— ^ = 0. (18) 

dx dx 

We introduce the following scaled dimensionless coordinate, 

x = E y. (19) 

The model is therefore characterized by the following parameters, 

n {x m k B T) 2 Ar 0+ r 0+ _ 

To= — D~E 2 — ' V = E' X = ^( a P+ + a P-)' P = -^( a P+~ a P-h ( 20 ) 

The coupled differential equations then read, 

da d 2 a dS . . 

-Ty + J? + X = l }i + <%' (21) 

dS ^ d 2 5 \§ + p ^° ^22) 
dy dy 2 dy : 

and a(y) and 5(y) are subject to the BC, 

,(, = ,0=0, d ^p2> =0 , % = Sl ) = 0, to*) =0 . (23) 
dy dy 

The solution of this system gives, 



3 



tf( y ) = 5] ai exp( ft y) + ^-M (24) 



i=i 



(A - /fy) 
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qi + <?2 + <?3 = 2, 

qm + q-iqs + qm = l - A - if, 
q^qs = -{X- f3rj). 

The a, are given as follows, 
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(25) 



a i 



qz{ql - A) - q 2 {qt - A) ) exp((g 3 + q 2 )y ) 



D S (X-(3 V ) 

~Jj~( X(y \_ ~ P) ( ?3 ex P(<?32/o + q2Vi) - q2 exp(q 3 y 1 + q 2 y )^j . (26) 

Here a 2 and a 3 are determined by cyclic permutation of the indices of the q^ from a 1; 
and D s is given by 

D s = gi<? 2 (<?2 ~ <?i)(<?3 - 1) exp((gi + q 2 )y + q 3 y{) + cycl.perm. (27) 

At y , we have: 

3 

°"(j/o) = -Y]ai(-l + ?i ) exp^o) - (ft ) exp(ft2/i) 

LV ?i y V qj 



(f3-X V ) + rXtf-Xr,) (3 



(yi - yo ) 



r](X — (3rj) l7] (X — (3r]) r]. 
Combining Eq. (I13|ll5p . the mean exit time for the Brownian particle follows: 
a(y )+5(y ) a(y ) - 5(y ) 



2a+ 



+ 



2a~ 



(2f 



(29) 
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Figure 6. (left) MFPT as function of a + and a . (right) current as function of a + 
and oT . 



The effect of shot noise on ion channel conduction is related to the arrival of ions at 
the channel's mouth. We calculated the channel MFPT as a function of the two flipping 
rates. Since the channel has a high affinity for the ions and therefore slow them down, 
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we have used a smaller ion diffusion coefficient inside the channel equal to 1.33 x 10~ 10 
m 2 s _1 . The results of the calculations are shown in the Figs. [5] and [61 It is clear from 
Fig. [5] that the particle takes longer to cross AEq as compared to AE±. There is a 
fast drop of the MFPT as the flipping rate a + increase. The current I = e/r is also 
presented as function of a + . A more general view of the dependance of the MFPT 
and the current on the two rates is shown in Fig. [61 We emphasize that the obtained 
transition rates are essentially non-equilibrium. In particular, current saturation effect 
can be observed as a + increases in a wide range of parameters. This is in accordance 
with experimental observations [33] of the current saturation at high concentrations. 

5. Conclusion 

In summary, we have introduced a Brownian dynamical model of ionic transitions 
through a channel, taking into account charge fluctuations at the channel mouth and the 
fluctuations of the channel walls. The statistical properties of the charge fluctuations 
are reconstructed from 3D Brownian dynamics simulation of multi-ion motion in the 
bulk solution. It is shown that distributions of ion arrival times at the channel mouth 
are exponential. It is further shown that these charge fluctuations strongly modulate 
the potential barrier for the conducting ion at the selectivity site due to amplification 
of electrostatic interactions in long narrow channels of low dielectric constant. These 
findings have allowed us to model the mean ion transition time through the channel 
as an ionic escape from the potential wall at the selectivity site induced by thermal 
fluctuations and by modulation of the height of the barrier by stochastic Poisson 
processes. The derived model is a Brownian dynamical model of the "knock-on" 
mechanism of the type [40l |4T]. Our model allows for analytic estimation of transition 
probabilities in the presence of charge fluctuations, i.e. it allows for analytic estimation 
of correlations between bulk concentrations and ion currents in charged narrow channels. 
In particular, it demonstrates the effect of current saturation due to ion concentrations 
in the bath. The model is essentially of a non-equilibrium nature. This last point is 
worth emphasizing because traditional approaches assume equilibrium rates of hopping 
between the sites. 

We note that the model takes into account the wall fluctuations. The latter feature 
is very important e.g. for an analysis of the tightly correlated motion at selectivity 
site of the type discovered in KcsA [42J. Our model allows for analytic estimation 
of transition rates in the presence of oscillations of channel walls, using our earlier 
results [43] on escape from periodically driven potentials using the method of logarithmic 
susceptibility [44] as will be discussed in details elsewhere. Such an escape process, 
assisted by the periodic modulation of the potential barrier by the wall oscillations, 
can result on its own in selectivity between alike ions due to the difference in their 
diffusion coefficients [431 14*6]. However, ultimately the selectivity of the channel has to 
be incorporated into the model by taking into account the effects of hydration [37] (see 
also [20] for a discussion of how hydration effect can be incorporated in our model). It 
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is also worth mentioning that the model that takes into account fluctuations of the wall 
may account for the dissipation of energy in the channel, for self-induced acceleration 
of the transition rate of the ion through the channel, and for coupling of the ion motion 
to the channel gating mechanics. Indeed, in this model, part of the energy induced 
by a very strong Coulomb interaction between the charged site at the channel wall 
and the moving ion is stored as the energy of vibrational modes. The later energy 
is only partially dissipated by the protein phonon modes. The remaining energy can 
now be used to modulate the potential barrier for the moving ion in a self-consistent 
manner to accelerate its transition through the channel. It can also be used to assist 
the conformational changes leading to the channel gating. 

The work in progress contains a plethora of unsolved problems. The immediate 
extension of the proposed model will be to include more then two levels for the potential 
at the selectivity site, taking into account positive, neutral, and negative charge at the 
channel mouths at both ends of the channel. The model can be further refined by 
including the estimates of the return times corresponding to a return of the ion in the 
channel to the initial bulk solution. 

Following the discussion above we can also formulate the following general unsolved 
problems in ion channels. We believe that they can be tackled by extending the model 
introduced in this paper: 

(i) The role of the membrane fluctuations. 

(ii) The role of the hydration potential. 

(iii) The role of additional binding sites outside the selectivity filter. 

(iv) The energetics of the ion transition including energy relaxation due to the coupling 
to the protein phonon modes (wall oscillations). 

(v) The coupling of the ion-wall interaction to the gating mechanism. 

In all of these, noise and dynamical effects seem to play a crucial role that is only starting 
to be elucidated. Our preliminary research shows that the model can be extended to 
take all of these effects into account. 
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